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We investigate how large baryon densities (and possibly high temperatures) may induce spon- 
taneous parity violation in the meson sector of QCD. The analysis at intermediate energy scales 
is done by using an extended er-model lagrangian that includes two scalar and two pseudoscalar 
multiplets and fulfills low-energy QCD constraints. We elaborate on a novel mechanism of parity 
breaking previously proposed by the authors based on the interplay between lightest and heavier 
meson condensates, which therefore cannot be realized in the simplest o model. We emphasize that 
the mechanism proposed here differs from the old idea of pion condensation advocated originally by 
Migdal. The results are relevant for an idealized homogeneous and infinite nuclear (quark) matter 
I /* ! \ where the influence of density can be examined with the help of a constant chemical potential. The 

model is able to describe satisfactorily the first-order phase transition to stable nuclear matter, and 
predicts a second-order phase transition to a state where parity is spontaneously broken. We argue 
that the parity breaking phenomenon is quite generic when a large enough chemical potential is 
present. Current quark masses are explicitly taken into account in this work and shown not to 
change the general conclusions. We expect that our approach will be adequate for dense nuclear 
matter of a few normal densities where quark percolation does not yet play a significant role. 
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. I. INTRODUCTION 

vo 

■ Emergent parity violation for sufficiently large values of the baryon chemical potential (and/or temper- 

ature) has been attracting much interest during several decades (see reviews [l[). It is actively looked for 
in both in dense hadron matter (in heavy ion collisions at intermediate energies and neutron stars) and in 
strongly interacting quark-gluon matter ( "quark-gluon plasma" in heavy ion collisions at very high energies 
and in quark stars). At finite baryon density pion condensation is a plausible possibility (conjectured in 
nuclear physics long ago in [2j. Although an attractive idea, the consensus reached by workers in the field 
is that it is not possible in simple models of pion-nucleus interaction. However in this work this idea will be 
vindicated and in fact claimed to be a very realistic possibility in more elaborated models. 
5^ 1 More recently the phenomenon of parity breaking was assumed to be present in meta-stable nuclear bubbles 

with non-zero axial charge created in hot nuclear matter [3( and/or in the presence of a strong background 
magnetic fields [3, Q . It was also shown Q that the associated axial chemical potential causes the distortion 
of energy spectrum of photons and vector particles (p and lo mesons) due to Chern-Simons interaction and 
as a consequence of parity breaking a remarkable mass splitting of their left- and right-handed polarizations 
appears. This phenomenon is however theoretically somewhat different in its origin and unrelated to the 
present paper, which is devoted to search for phase transitions with generation of neutral pseudoscalar 
isospin-one condensate. 

Parity violation may in fact accompany the transitions to open color phases like the color- flavor locking 
(CFL) [tJ or superconducting ones 8|, but these regions are beyond the scope of the present paper. 

Some time ago it was proved in Q that parity and vector flavor symmetry could not undergo spontaneous 
symmetry breaking in a vector like theory such as QCD at normal vacuum conditions at zero chemical 
potential 10]. Finite baryon density however may result in a breaking of parity invariance by simply 
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circumventing the hypothesis of the theorem. Indeed the presence of a finite chemical potential leads to the 
appearance of a constant imaginary zeroth-component of a vector field and the conditions under which the 
results of [9( were proven are not fulfilled anymore. 

Can this possibility of parity breaking be realized in nature and what would then be its experimental 
manifestation? In order to answer the first question we might appeal to lattice QCD for help and in fact 
this possibility has been studied intensively for quite some time. However, finite density simulations are 
notoriously difficult (the fermion determinant is in general complex for a non-zero chemical potential ^i) and 
one has to resort to more involved techniques such as determining the phase of the determinant separately, 
Taylor expansions in /i or analytic continuation to imaginary chemical potential. On top of that Pauli 
blocking is at work for large values of fj,, etc. For a comprehensive review of the ongoing issues, see for 
instance fl3l - [l5j . Thus the lattice results for sufficiently large values of the baryon chemical potential (where 
the effect is expected to appear) are not known rigorously yet [16| . 

In this work we shall attempt to explore the interesting issue of parity breaking employing effective 
lagrangian techniques, useful to explore the range of nuclear densities where the hadron phase still persists 
and quark percolation does not occur yet. Our effective lagrangian is a realization of the generalized linear a 
model, but including the two lowest lying resonances in each channel, those that are expected to play a role 
in this issue. This seems to be the minimal model where the interesting possibility of parity breaking can 
be realized. The use of effective lagrangians is also crucial to answer the second question of interest, namely 
how would parity breaking originating from a finite baryon density eventually reflect in hadronic physics. 
This work is an extension of the preliminary results presented in [l8l . [l9| . 

An oldest, pre-QCD attempt to give the lagrangian description for two multiplets of scalar and pseudoscalar 
mesons was done in [20( with a reduced set of operators and a chiral symmetry breaking (CSB) pattern not 
quite compatible with QCD. In turn, we have been rather inspired by our previous works on extended 
quark models (2l| - [23j where two different schemes with linear and non-linear realization of chiral symmetry 
were adopted to incorporate heavy pions and scalar mesons within an effective quark model with quark 
self-interactions. A certain resemblance can be also found with the model [24| where two SU (3) p multiplets 
have been associated with two-quark and four-quark meson states although we do not share the assumption 
in [24j concerning the dominance of four-quark component in radially excited mesons. The model in (2|| is 
also of relevance in studying of vacuum for extended a models. 

The paper is organized as follows. In section 2 the bosonization of QCD quark currents in the color- 
singlet sector is discussed and the ingredients of the generalized a model are indicated. Some subtleties 
associated to the choice of the low-energy effective hadronic theory are identified. Too simple models are 
not rich enough to explore all the different phases that the presence of parity violation opens for us. In 
section 3 we introduce the a model with two multiplets of isosinglet scalar and isotriplet pseudoscalar fields. 
The effective potential for two multiplets of scalar and pseudoscalar mesons is obtained and the mass-gap 
equations and second variations at the minima are derived. We impose on the model the conditions for it not 
to lead to spontaneous breaking of parity for vanishing baryon chemical potential. This condition bounds 
a combination of several coupling constants. In section 4 the freedom of reparameterization is exploited to 
reduce the number of relevant constants so that the mass-gap equations allow leads to physically inequivalent 
solutions. 

In section [V] we shall introduce the finite chemical potential and temperature and see how they modify 
the effective theory and the vacuum state. Temperature and baryon chemical potential appear through the 
one-quark loop free energy; the only one that needs consideration in the large JV C limit as the contribution 
of meson loops is neglected in the mean field approach. Different thermodynamic properties of the model 
are obtained for non-zero temperature and chemical potential (pressure, energy density and entropy). 

In the subsequent sections we ignore temperature effects and concentrate on the finite density case. In order 
to describe adequately the saturation point transition to stable nuclear matter we supplement the effective 
lagrangian with a u meson coupling to the isosinglet quark current. It is known for some time (26j that ui 
mesons contribute to the repulsive part of nuclear forces and thereby they support the formation of stable 
nuclear matter. Indeed, the correct first-order phase transition at the saturation point is found and we prove 
that the chiral collapse affecting the simplest a models and/or the one-mutiplet Nambu-Jona-Lasinio (NJL) 
models does not occur here. Once the saturation point is established one can define (in) compressibilities 
and perform the matching between quarks and nuclear matter. 

In section 6 we investigate the emergence of spontaneous parity breaking (SPB) phase. The mass gap 
equations and critical lines for the parity breaking phase transition are derived and corrections beyond the 
chiral limit taken into account to the leading order in quark masses. An interesting question examined is on 
how the saturation point meets spontaneous parity breaking. 

In section 7 we study the approach to the SPB phase transition both in the chiral limit and beyond. It 
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is established that this phase transition is of second order. The second variation of the effective potential 
leads to the mass spectrum for scalar/pseudoscalar states in this phase. We remind that in the SPB phase 
strictly speaking there are no genuine scalar or pseudoscalar states as each of massive states can equally 
well decay into two and three scalars. In section 8 the kinetic terms are considered in order to determine 
the physical masses of mesons. In the SPB phase the masses of the lightest states are calculated beyond the 
chiral limit and the existence of two massless Goldstone bosons verified. We summarize our findings and 
propose possible signals of the SPB phase in the conclusions section. 

The range of intermediate nuclear densities where our effective lagrangian is applicable is of high interest 
as they may be reached in both compact stars [13] and heavy-ion collisions 28]. Its relevance can be 
qualitatively motivated by the fact that at substantially larger densities typical distances between baryons 
are shrinking considerably and meson excitations with Compton wave lengths much shorter than a pion one 
start playing an important role. 

There are previous studies dealing with the problem of strong interactions at zero temperature and finite 
chemical potential: depending on a value of nuclear density, a variety of methods are involved from using 
meson- nucleon [EHil or quark -meson 27, 30] lagrangians for low-dense nuclear matter to models of the NJL 
tvpe(3l| for high-dense quark matterQ^At very large chemical potentials, when diquark condensation gives 
rise to a new CFL phase one can guess Q an effective lagrangian describing the lightest degrees of freedom 
and the conditions for SPB due to pion or kaon condensation. 

Recent papers are more concerned with case of non-zero chemical isospin potential f32| (where comparison 
with lattice results is possible). In our view these works are somewhat incomplete and hence they may not 
lead to conclusive answers for intermediate nuclear densities when quark percolation does not yet occur but 
pion-nucleon theory is certainly insufficient to give an adequate description. Although the issue of SPB in 
hadronic phase has been touched upon in the pion-nucleon theory[l|, [2^ . l33j and in NJL models [34| the 
validity of the models used is not quite clear for intermediate nuclear densities. The reason is discussed in 
the next section: they are not rich enough to explore the subtle phenomenology involved. There are also 
promising recent attempts [351 ] to describe the large baryon density influence on hadron properties. 



II. BOSONIZATION OF QCD IN THE COLOR-SINGLET SECTOR 



In order to elaborate an effective lagrangian for meson states starting from QCD we revisit the properties 
of color-singlet (quasi)local quark currents in the QCD vacuum with spontaneously broken chiral symmetry. 
This phenomenon emerges due to a non-zero value of quark condensate (qq) and can be associated to the CSB 
scale Aqsb ~ lGeV(= A for brevity). This CSB due to quark condensation makes the quark bilinears to 
be interpolating operators for meson fields (in the limit of large colors). In particular, a scalar/pseudoscalar 
quark density effectively describes the creation or annihilation of scalar/pseudoscalar mesons 

qq(x) ~ A 2 J2 Z^ix)- q l5 r l q{x) ~ A 2 £ zj 1) Ilj(x), (1) 
I I 

with the normalized meson fields , IT| describing the families of resonances with the same quantum numbers 
but increasing masses (radial Regge trajectories) [36] and the set of normalization constants Z\ X \ The 
constituent quark fields are denoted as q, q and t 1 stand for Pauli matrices. At this stage we consider the 
chiral limit of zero current quark masses. Accordingly the global chiral covariance of quark operators (JT|) is 
transmitted to the set of boson operators leading to an equal normalization of scalar and pseudoscalar fields. 
This is a basic framework of linear sigma models [37| (see next section). In this paper we restrict ourselves 
with consideration of two light flavors related to u, d quarks and therefore the approximate chiral symmetry 
of the quark sector is SU{2) L x SU(2) R . 

Keeping in mind confinement we have retained in ([1} only the one-resonance states as leading ones while 
being aware of that the total saturation of quark currents includes, of course, also multi-resonance states. 
Thus we use the large- N c approach where resonances behave like true elementary particles with zero widths 
and multi-resonant states can be neglected. 

Let us comment a bit more on the previous relation. On the left-hand side one sees an operator of canonical 
dimension 3 whereas on the right-hand side one finds field operators of canonical dimension 1 . This drastic 
change in dimensions is a consequence of CSB and it modifies the dimensional analysis of what must be 
included into an effective lagrangian. More exactly, in order to replace the non-perturbative regime of QCD 
at low and intermediate energies by a hadron effective lagrangian one has to apply this dimensional counting 
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in the CSB phase [22j to all possible combinations of color-singlet operators arising in the chiral expansion 
in inverse powers of the CSB scale A. 

To be specific the chiral invariant local operators playing the leading role in the low-energy effective 
lagrangian for meson dynamics are 



^2 - qj 5 nqq~f 5 T l q] ~ A 2 £ 4nP' S ™ + n?IL., m ] + 0(1); 



1,111 



^[{qq? - q^nqq^q} 2 ~ £ Z% nr Wm + HfHi.m] [(E„Er + 11^11^] + O(^); 

Z,m,n,r 

^(rf - ^(^5)^75^)] ~ X! ^ } [^S;^s m + d^u^} + O(^), (2) 

where the matrices z{f,^„ r ,2^ must be symmetric under transposition of indices in order to provide 
global chiral invariance. The terms quadratic in scalar fields trigger (if their coefficient is of the right sign 
an instability that leads to CSB in the effective meson theory (see. e.g. [HI). 

The above set of operators is not complete and can be extended with the help of form factors that are 
polynomials in derivatives [22J. For example, using the same CSB scale A one can add into the effective 
quark lagrangian the vertices built of the elements 

^^W^T^'^'W; ^^(^A 2 ^ 1 ^*), (3) 
/ i 

which may give numerically comparable contributions for several k [22j . 



III. GENERALIZED SIGMA-MODEL 
A. Effective potential for two multiplets 

The simplest hadronic effective theory is the linear sigma-model of Gell-Mann and Levy [37J, which contains 
a multiplet of the lightest scalar a and pseudoscalar ir a fields. Spontaneous CSB emerges due to a non-zero 
value for (er) ~ (qq)/F 2 . Current algebra techniques indicate that in order to relate this model to QCD 
one has to choose a real condensate for the scalar density, with its sign opposite to current quark masses, 
and avoid any parity breaking due to a v.e.v. of the pseudoscalar density. The introduction of a chemical 
potential does not change the phase of the condensate and therefore does not generate any parity breaking. 
This is just fine because in normal conditions parity breaking does not take place in QCD. However, if two 
different scalar fields condense with a relative phase between the two v.e.v. 's the opportunity of spontaneous 
parity breaking may arise. 

Let us consider a model with two multiplets of scalar ( Bj) and pseudoscalar ( 7r° ) fields 

H^ajl + iTTj, i = l,2; H^) = (5* + (tt]) 2 )I, (4) 

with TTj = TTjT a with r a being a set of Pauli matrices. We shall deal with a scalar system globally symmetric 
respect to SU(2)l x SU(2)r rotations in the exact chiral limit and next consider the soft breaking of 
chiral symmetry by current quark masses. We should think of these two chiral multiplets as representing 
the two lowest-lying radial states for a given J PC . Of course one could add more multiplets, representing 
higher radial and spin excitations, to obtain a better description of QCD, but the present model, without 
being completely realistic, already possesses all the necessary ingredients to study SPB. Inclusion of higher- 
mass states would be required at substantially larger densities when typical distances between baryons are 
shrinking considerably and meson excitations with Compton wave lengths much shorter than a pion one 
start playing an important role. 

Let us define the effective potential of this generalized sigma-model. First we write the most general 
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Hermitian potential at zero /i, 



\ tT | - 12 H j A J* H k + XiiHlHt) 2 + \ 2 {HlH 2 ) 2 + \zH\HxHIH2 

+ h 4 (HlH 2 HlH 2 + HlHxHiHx) + ^{H\h 2 + H^H^ 

+ \\ & {H\H 2 + fftffOfftfla j + O(^f), (5) 

which contains 9 real constants Ajfe, A^; A — 1, ... ,6 . However this set of constants can be reduced (see 
sect II V[) . QCD bosonization rules in the large N c limit prescribe Ajk ~ A^ ~ N c . The neglected terms will 
be suppressed by inverse power of the CSB scale A ~ 1 GeV. If we assume the v.e.v. of Hj to be of the 
order of the constituent mass 0.2 -j- 0.3 GeV, it is reasonable to neglect these terms. 

One could add five more terms (breaking parity manifestly): an imaginary part of A12 with an operator 

itr{(ff^ff 2 -#+#0}, (6) 

and three more operators 

itr I H\ H 2 h\ H 2 - H\ H x h\H-X, itr j [h\ H 2 - h\ H x ) h\ H x X , 

itx[{H\H 2 - HlHx)HlH 3 } . (7) 

As well as an operator with disconnected trace seems to complete the full set of operators 

(ir[H\H 2 )) 2 +(ir[H 1 Hl)) 2 . (8) 

However after specifying the v.e.v. (Hi) = (<ti) one can use the global invariance of the model to factor out 
the Goldstone boson fields with the help of the chiral parameterization 



Hx(x) = cr 1 (x)£, 2 (x) = ci(x)exp (^~J^J 5 H 2 {x) = £(x)(a 2 (x) + in 2 (x)j £(x) , 



(9) 



which differs from eq. (U) in notation. This kind of parameterization preserves the parities of u 2 (x) and 
n 2 to be even and odd respectively in the absence of SPB. Then the contribution of the four additional 
operators © , vanishes identically, whereas the operator (|5]) turns out to be a combination of operators 
with constants A3, 4. Finally the potential (|5|) is further simplified to 

2 2 
V eS = - 12 ^A^ fe -A 22 (7rf) 2 + A 2 ((7r^ 2 ) + K) 2 ((A 3 - A 4 )a 2 + X 6 aia 2 + 2A 2 a 2 ) 

j,k=i 

+\icjI + \ 2 o\ + (A 3 + \±)o\a\ + X 5 afa 2 + Ae^of. (10) 

The current quark mass m q corresponds to the average of the external scalar sources Mj(x) = Sj[x) + 
iT a Pj(x), namely, (Mj(x)) — — ^djin q and thus the relevant new terms beyond the chiral limit can be 
produced with the help of the formal replacement Hj — > Cjm q in all quadratic and quartic operators included 
in ([5]) and by adding these new terms with new constants into the effective potential. We will consider the two 
flavor case and retain the softest terms linear in Hj and m q , thereby neglecting terms cubic in scalar fields 
exploiting the non-linear equivalence transformation Hj — > J^k 1 m=i 2 bjkimHkH[ H m . With an appropriate 
choice of external scalar sources only one type of Hermitian structures is possible, Ylj=i 2 ^ r (MjHj + HjMj) 
, i.e. we add two new terms to our effective potential (|10p 



1 

-m q tr 



d l {H l +Hl) + d 2 {H 2 + Hl)\. (11) 
Making use of our chiral parametrization of the fields Hj through the chiral field 



9 |7T?| T a 7T? |7T?| , s 

U^e = cosL^ +l — fsin^i, (12) 



() 



one derives the following extension of the effective potential (fTU)) 



&V eS (m q ) = 2r 



-(did + d 2 a 2 ) cos i-fi + (fe-p-rf sin i-fi 

^0 Fl I ^0 



(13) 



The effective potential (fTU|) . (ITUl) will be used to search for CSB and for the derivation of meson masses. 

B. Mass-gap equations and second variations 

Let us now investigate the possible appearance of a non-zero v.e.v.'s of pseudoscalar fields. In order not 
to violate charge conservation, we may expect in the pseudoscalar channel the condensation of the neutral 
isospin components only, 

(tt?) = (n°)S 0a , (4) = p6 0a . (14) 

Now one obtains four mass-gap equations as, unlike in the chiral limit, the pion field becomes in principle 
observable 

(tt ) 

2(AuCTi + Ai 2 cr 2 ) + 1m q d x cos ^— - = A\\o\ + 3X 5 a 2 a 2 

+2(A 3 + \ 4 )ai<4 + \ &a \ + p 2 / 2 ^ 3 _ x ^ ai + A \ ^ ( 15 ) 

(tt ) 

2(Ai 2 cri + A 22 cr 2 ) + 2m q d 2 cos ^— - = \ b a\ + 2(A 3 + \±)o\o 2 

+3\ 6 <7ial + A\ 2 a\ + p 2 (\ e ai + 4A 2 o- 2 J , (16) 
(tt ) (tt ) 

(di<7i + d 2 o- 2 ) sin — — = -d 2 pcos — — , (17) 
-ro -ro 

. (tt ) m q dl (vr ) 

~m q d 2 sin — — = p j- r cos — — 

F Q (ditTi + d 2 a 2 ) F 

= p(- A 22 + (A 3 - X 4 )cr 2 + \a x a 2 + 2\ 2 a\ + 2\ 2 p 2 ) , (18) 



As follows from eq. (fl7|) , the equality p — entails (7r°) = if dia% + d 2 a 2 ^ and d 2 ^ 0. However, as will 
be seen below, the combination di<xi + d 2 a 2 is related to the quark condensate 

d\G\ + d 2 a 2 = -(qq) > 0, (19) 

hence, this combination cannot be zero. For d 2 = one has always (tt ) = and the parity breaking 
pattern remains the same as for the massless case. For d 2 ^ both pseudoscalar v.e.v. (tt ) and p can 
arise simultaneously only. To avoid spontaneous parity breaking in then normal vacuum of QCD, it is thus 
necessary and sufficient to impose 

771 d 2 

(A 3 - Xi)al + \g x g 2 + 2\ 2 a\ > A 22 + j-, (20) 

Since QCD in normal conditions does not lead to parity breaking, the low-energy model must necessarily 
fulfill (|20|. 

In the usual vacuum state the necessary condition to have a minimum for non-zero Oj (for vanishing 
p), equivalent to the condition of having CSB in QCD, can be derived from the condition to get a local 
maximum (or at least a saddle point) for zero aj . At this point the extremum is characterized by the matrix 
— Ajk in ([5]). It must have at least one negative eigenvalue. This happens either for DetA > 0, tr{A} > 
(maximum at the origin) or for DetA < (saddle point at the origin) . The sufficient conditions follow from 
the positivity of the second variation for a non-trivial solution of the two first equations l|15p. ([To]) at p = 0. 
The matrix containing the second variations for the scalar sector is 

-Vg )a = -An + 6Aia? + 3A5CJK72 + (A3 + M)<J 2 2 , 

V} 2 ^ = -2A 12 + 3A 5 ct? + 4(A 3 + \ 4 )<Ji<J2 + 3A 6 af , 

^V 2 (2)a = -A 22 + (A 3 + A 4 )a? + 3A 6( j 1( 72 + 6A2CT2 2 . (21) 



7 



In turn, the nonzero elements of the second variations in the pseudoscalar sector are 

1 (V^\ ab - S ab m dl(Jl + d2 ° 2 

-^v n j —o m q , 



^(V£ )n ) ab = S ab (- A 22 + (As - + X 6 a lf y 2 + 2X 2 a(] . (22) 



F ' 

The required conditions are given by tr > and DetV"( 2 ) > 0. For positive matrices it means that 

^ 2)ff >0; V^>0. (23) 

The eigenvalues of (|22|) eventually give the masses squared of tt, n' mesons and thereby must be positive 
according to the inequality (|2T))) . The latter corresponds to the positivity of the determinant. From the 

positivity of matrix element it follows that the combination d\(j\ + d 2 a 2 > 0. In fact, to the leading 

order in m q the masses of tt and n' mesons are proportional to Vti and , respectively, whereas the 
off-diagonal matrix element V$' K contributes to the second order in m q . The requirement to have a positive 
determinant of the matrix vj£ is exactly equivalent to (l20l) . 

The two set of conditions, namely those presented in eq. (|20l) and in eq. (1231) represent restrictions that 
the symmetry breaking pattern of QCD imposes on its low-energy effective realization at vanishing chemical 
potential. 

One can easily find the correction linear in m q to the vacuum solution in the chiral limit 

a j (m q ) ~ a, (0) + 2m g AJ ; (24) 



Using these equations the corrections to the masses of scalar and heavy pseudoscalar mesons can be derived 
straightforwardly. In particular, for scalar mesons the corrections to the mass matrix are 



= 2m q £ d jV ^(V^):y m (25) 

j,m=X,2 

whereas in the pseudoscalar sector 

AV^ = 2m q £ ^4 2)7r (F^)7^ m . (26) 
j',m=l,2 



The latter term saturates the current quark mass correction for heavy pseudoscalar meson masses. 



IV. REDUCTION OF COUPLING CONSTANTS AND EXTREMA OF U eff 

Let us investigate how many extrema the effective potential possesses for different values of the coupling 
constants. In this section we take the chiral limit m q = for simplicity. It turns out that when the chemical 
potential and temperature are zero one can eliminate one of the constant in the effective potential by a 
redefinition of the fields. Indeed, one can change the variable 

H 2 = aHi + 2 , (27) 

using a linear transformation with real coefficients a, ft (to preserve reality of <Tj,7r°). With the help of 
this redefinition one can diagonalize the quadratic part in ((3J) and make its coefficients equal An = A 22 = 
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det A/A 22 = A. Then 

2 

]T tr {i?]A jfc H fe } = Atr {fl^ + iT^} . (28) 

3,k=l 

A further reduction of the coupling constants affects the dependence of free energy on finite chemical potential 
and temperature (see below), but it can be implemented when both external control parameters vanish; 
namely we perform the following orthogonal rotation of two fields 

Hi = cos c/)Hi + sin (f>H2, H2 = — sin cj)Hi + cos 4>H^. (29) 
Then the coefficient in the operator (H[H2 + h\Hi)h\H\ becomes equal to 

X5 . , = A 5 -2(A 3 + A4-2Ai)tan0-3(A 5 - A 6 )tan 2 </> 
cos 4 cp 

+2(A 3 + A 4 - 2A 2 ) tan 3 <j) - A 6 tan 4 <f> = V\ 5 (t&n(f)). (30) 

One can always fix A§ < by reflection of H2 ■ Then if A5 < then V\ 5 (0) < but evidently for tan<^> ^> 1, 
V\ s (t&n<t>) ~ — Agtan 4 (/) > and therefore the equation 7- > A 5 (tan0) = has at least one (positive) real root 
. In the complementary region A5 > and therefore P\ 5 {0) > 0. In this case one can look at tan</> = ±1 
where 

Vx 5 (±1) = -2(A S - A 6 ) ± 4(Ai - A 2 ), (31) 

so that one of these combinations is negative. Again the comparison with the asymptotics allows to conclude 
that there is a real root for V\ s (teLTi(j)) = 0. Thus for any sign of A5 it can be eliminated by a proper rotation 
of scalar fields. 

Let us take the basis of operators with A5 = . Then, after renaming the fields 

V eS = -A(Vx) 2 + M 2 ) + A 2 ((7r 2 a ) 2 ) 2 + (7r 2 a ) 2 ( - A + (A 3 - A 4 )a 2 + A 6 a lCT2 + 2A 2 a 2 ) 

+A10- 4 + A2CT2 + (A 3 + \±)o\o\ + Ka x a\. (32) 

This potential simplifies the mass gap equations and second variations in order to investigate their solutions 
analytically. The effective potential must provide the familiar CSB at normal conditions (p = T = 0). Thus 
in the chiral limit there are at least two minima related by the symmetry rotation 01.2 — > —01,2 and one 
maximum at the origin. This is implemented by assigning a real singlet v.e.v. (a\) > to H\ thereby 
selecting one of the minima. 

In this section we shall assume A5 = in order to determine the different vacua of the theory at zero 
temperature and chemical potential. 



A. Search for the extrema of the effective potential 



In the parity symmetric case the second (|16l) reads 



3 

cr 2 (-A + (A 3 + A 4 )f7 2 + -Ka X 02 + 2A 2 <r|) = 0. (33) 



One solution is 



4 0) =0, (ai 0) ) 2 ^, (34) 

from eq. (fT5|) . For stable solutions Ai > and therefore A > 0. 

Another set of solutions cr^; a — 1,2,3 comes from eq. (f3"3"| for er 2 ^ . With a combination of the 
mass-gap eqs. (|15|) and (|33j) one can decouple the equation in terms of the ratio t = 172/ci, 

V 3 (t) = t 3 -at 2 -bt + c = 0, 

2((A 3 + A 4 )-2A 2 ) 2((A 3 + A 4 )-2A0 
a = ^ 6 = 3, c=^ '-, (35) 

— Afi — Ak 
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where the sign is fixed for A§ < and c > as it is shown in the next subsection. As the order of the 
equation is odd there may one or three real solutions. Because 7^3 (0) > 0, 7-3(0) < one concludes that one 
of the solutions is negative. Let us analyze the extrema of V$(t) 



Vs(t)=0 — >t 2 --at-b = 0, t± = -a±J-a 2 + b,t + >0,t-<0. (36) 
3 3 V 9 

All together it means that a negative solution tW < always exists and (if any) two more solutions are 
positive, < 2a 3 ), and separated by a minimum. Therefore the existence of two positive solutions is 
regulated by the sign of Vs(t+). They exist if 

Ps(t+) _ -A,,* - ^ + e . i^^l = <=-§(f + ^) 1 < 0. (37) 
Evidently it takes place for positive a, i.e. for 

(A 3 + A 4 ) > 2A 2 . (38) 
Later on we will see that in order to implement a first-order phase transition to stable nuclear matter we 




FIG. 1: Extrema of effective potential in the reduction basis: the maximum is placed in the square, up to four 
minima are located in the circles and the corresponding four saddle points are depicted by the lentils. The existence 
of two solutions t^ 2 \t^ 3 ' with positive values of Oj is governed by condition i|37)) . Which one corresponds to the true 
minimum depends on the actual value of the phenomenological constants. 

need two minima for positive u\ which entails two more minima for negative o~j due to symmetry under 
o~j — > —crj and therefore four saddle points must connect them. Thereby, in the half plane of positive <t\ one 
has to reveal four solutions, namely, one is = and three for er 2 ^ so the condition (|3T[) should be 
satisfied in order to be able to describe the saturation point, see Fig.l . 

Let us recall that all the inequalities obtained in this section are referred to the field basis with fully 
diagonal Ay = ASij and with A5 = 0. 
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B. Selection of the minima 

In all cases the condition of minimum come from the positive definiteness of the matrix of second variations 
1. 



<Tl<7 2 
2 7T2^2 



For = they read 



-A + 6A 1( t 2 + (A 3 + A 4 )(t| >0, 
4(A 3 + A 4 )cri(T2 + 3A 6 cr 2 , 
-A + (A 3 + A 4 )cr 2 + 3A 6 ericr 2 + 6A 2 cr| > 0, 
-A + (A 3 - A 4 )cr 2 + A 6 cricr 2 + 2A 2 cr 2 > 0. 

(A 3 ±A 4 ) >2Ai,— > A 3 > |A 4 | 



(39) 



(40) 



for A > as it is required by the absence of chiral collapse (and the spectrum at the SPB point, see below). 
It gives support to the condition c > in the previous subsection. 

For cr 2 7^ one obtains a number of bounds on the solution from the second variation 



1 



Vg> n = (oi) s 



4Ai - l -\ & t z 



> 0, 



V™ = (<?i? 4(A 3 + A 4 )t + 3A 6 t 2 



~V {2) = (dl) 2 



2 7T27T2 



-Aat + 4A a f 
— 2A 4 — 2^6* 



>0, 
> 0. 



(41) 



Evidently if A 4 > then <j<x > — > t > 0. The remaining bound must come from the positivity, det V^ 2 ' > 



FINITE TEMPERATURE AND BARYON CHEMICAL POTENTIAL 



A. Coupling the effective lagrangian to the environment 



We are building a model of meson medium starting from the quark sector of QCD. Its thermodynamical 
properties and relationship to a dense baryon matter will be examined with the help of thermodynamical 
potentials derived from the constituent quark model in the mean field (large N c ) approach. This gives a 
prescription to connect the properties of quark and nuclear matter and estimate the parameters of our model 
to reproduce meson phenomenology and the bulk characteristics of nuclear matter such as binding energy, 
normal nuclear density and (in) compressibility. 

The meson degrees of freedom present in our model appear after bosonization of QCD in the vacuum. The 
effects of infinite homogeneous baryon matter on the effective meson lagrangian is described by the baryon 
chemical potential ^i, which is transmitted to the meson lagrangian via a local quark-meson coupling (in the 
leading order of chiral expansion yit 2 /A 2 ). In turn, in the large N c limit one can neglect the temperature 
dependence due to meson collisions. The temperature T is induced with the help of the imaginary time 
Matsubara formalism for quark Green functions [39| 



(2n+ 1)tt 



1 

kf' 



(42) 



For real physics with 3 colors this approximation to thermal properties of mesons is expected to be less 
precise as meson loops contribute substantially to the thermodynamic characteristics. Nevertheless it should 
be sufficient to describe qualitatively the interplay between baryon density and temperature at the phase 
transition. 
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Without loss of generality we can regard one of the collective fields Hj, namely, Hi as the one having 
local coupling to quarks: this actually defines the chiral multiplet H%. The set of coupling constants in ((5|) 
is sufficient to support this choice as well as to fix the Yukawa coupling constant to unity. 

We notice that this specification of collective fields is compatible with the transformation (f27|) and therefore 
one can deal with the diagonal quadratic part of the potential. However the additional linear transforma- 
tion (I29[) would split the constituent mass in the quark Yukawa vertex into two fields qRHiq^ + h. c. — > 
qn(cos (pHi+sin </>i/2)<ZL+h. c. It means that a possible change of the basis used in sec. 4 to eliminate the con- 
stant A5 affects the chemical potential driver AV e ff(crx, /3, fx) —> AT4//(a/(cos cj> a\ + sine/) a 2 ) 2 + p 2 , /?,//). 
Thereby all the mass gap equations (|15|) - (I18I) obtain new contributions depending on T and p. 

To keep the consideration more tractable we do not reduce the coupling A5 — > but rather select the basis 
in which finite density and temperature is transmitted to the boson sector by means of 

AC = -(q R H iqL + q L H\q R ) — > -Q<nQ, (43) 

where Ql = Qr — tfqii stand for constituent quarks [23| . However, the results derived in the previous 
section on the different vacua for vanishing temperature and chemical potential remain obviously valid. 

Following this recipe the quark temperature and chemical potential dependence can be derived for mass 
gap equations - the conditions for a minimum of the effective potential. Namely, taking into account the 
transformation (f2T)f the first equation (fT5|) is modified to 

2Ac7i + 2m q di cos ^—^ = 4Ai<7? + 3X 5 a 2 a 2 + 2(A 3 -I- A 4 )<ria| + A 6 crf 
+p 2 (2(X 3 -X i )a 1 +X 6 a 2 ^ + 2Na 1 A{a 1 , p,j3), (44) 

with notation M = N , aN i f . In turn 

47T 

1 ^ r d?p 1 

A{crt,fi,P) = -Vexp(zw„77) / ■ — ^ =2 -[T = 0, /i = 0] 

P „ J T + p) z - E z 

d 3 p 



2ttE 



(f(E-p)-f(-E-p) + l) 



2 f dE ^E 2 - a 2 (f(E - M ) - f(—E - p) + 1) 



-2/ dEjE 2 -2 cosh(/3//) + exp(-/3£) 



where E = ^Jp 2 + erf and the limit 77 —> is to be taken after the Matsubara summation and the Fermi 
distribution function is introduced 

^ = c^MTT • (46) 

The corresponding one-quark loop contribution to V^g can be obtained by integration 

00 

AV e Q (a\,n,fi) = -2Af / dx xA(x, p, /3) 



--NJ dE(E -a,) cosh{M+cosHpE y (47) 



0-1 



so that it vanishes for very large a\ . 

Using the gap equations and (|4"4l , the value of the effective potential at its minima is given by the compact 
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expression 



m,, 



2~V 1 " ' r J 1 2""« 

1 



-(c!i(cri) + d 2 ((T2>) cos — h d 2 psuv 



Fq ' Fq 



AV eS «(Ti>,M„9) ~ -^(a!) 2 ^,/?)^^),/!^); (48) 
^r/ v2, / va^rf / N2 \ 1 /2cosh(/3/i)+cxp(-/3i : ;) 



AV eS - \m{ct 1 ) 2 A = -\m J dE(AE 2 — (ai) 2 )(^E 2 — (<ri) 2 



cosh(ftp) + cosh(ftE) 



P 



It should be emphasized that all the above results have corrections of O (/j, 2 /A 2 , cr 2 /A 2 ). In the remaining 
part of this section we simplify our analysis and take the chiral limit m q = 0. 

B. Thermodynamic properties of the model at T/0 

Thermodynamically the system is described by the pressure P, the energy density, e and the entropy 
density s. The pressure is determined by the potential density difference with and without the presence of 
chemical potential, dP — —dV 

(a 3 (p,ft),p,ft) = V eS (^(0,0))-^ ^(p,ft),p(p,ft),p,ft), (49) 

where the dependence of solutions of mass gap eqs. a 3 - and p on p and ft has been shown explicitly. The 
energy density is related to the pressure, density and entropy density by 

e = -P + N c pg B + Ts. (50) 

The chemical potential is defined as 

d BB e = N c p, (51) 

with the entropy and volume held fixed. The factor N c is introduced to relate the quark and baryon chemical 
potentials. Since e is independent of p 

oo 

d,P = N cQB , QB = -±8»V eS =^idE E (e 2 - a 2 ) 1 ' 2 .^^y (52) 
N c it 2 J \ J cosk(ftp) + cosh(ftE) 

ax 

In turn the entropy is defined as 

s = d T P = -d T V eS , Ts = ftd p V eS . (53) 
Let us relate the last two terms in (1481) and (l50l) 



Ay eff (a 1 , f i,ft)-^Afa 2 (p,ft)A(a 1 , f i,ft) = (1 - i^JAV^f (<n,lhP) 

= J(/i9 M - ftdp)AV eS (at,n,ft) = - ftd )V eS (a 1} p,ft) = -\{N c1 xq b + Ts), (54) 

The above equation allows to calculate the entropy in our model. 



C. Zero temperature and finite density 



In this subsection we consider the zero-temperature case and study the regime of chemical potentials 
comparable with the v.e.v. a%. At zero temperature T = the contribution from p to the effective potential 
is 



AV eff (oi,M) = y6»(/i-(Ti) 



2 > ■> , 2 1 ,:; 2 



0"! 



-a\ In : 



(55) 
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The total value of the effective potential at its minimum is 

A/' 



(56) 



Higher-order terms of the chiral expansion in 1/A 2 are not considered and the gluon condensate contribution 
is assumed to be hidden in the dynamical mass value. 
Accordingly in the first mass gap equation (j44| 



A(tn,n,P) =°°26([i-oi) / dE x /E 2 -af 



2 fl+ y/p 2 - a\ 

i - o-i in - 



0\ 



Then the second variation of effective potential is modified in the only element 



1 



^Vg )a = -A u + 6\ 1 ai+ 3A 5 a lCT2 + {X 3 + X 4 )o% 



VAf6{p - o x ) 



p,\J fj? — of — io l In 



2 , M + yV 



0\ 



The effective potential (I55l) . (|56p is normalized to reproduce the baryon density for quark matter 



Qb = -—df.AV^ (oi,/i) 



1 dV eS (») _ N 



- • f ^3 ly } ( .,2 _2( ,.\ 



N f 



N c dp 



3/2 



(57) 



(58) 



(59) 



where the quark Fermi momentum is pp = yj p 2 — of(p). 



VI. SPONTANEOUS PARITY BREAKING PHASE 

A. Mass gap and critical lines for the SPB transition 

Let us examine the possible existence of a critical point, in the chiral limit m q = for simplicity, where 
the strict inequality (1201) does not hold and instead for p > p cr u 

(A 3 - A 4 )cr 2 + XqOxOi + 2A 2 (<4 + p 2 ) = A. (60) 

After substituting A from (|60|) into the second eq. (fT5|) one finds that 

A 5 cr 2 + 4A 4 c7iCT2 + A 6 (of + p 2 ) = 0, (61) 



where we have taken into account that o\ 7^ 0. This, together with (|60|) completely fixes the v.e.v.'s of the 
scalar fields 0%^. If A2 = and/or Xq = equations tfSTTj) or (jfTLj) unambiguously determine the relation 
between o\ and <t 2 . Otherwise if A2A6 7^ these two equations still allow to get rid of the v.e.v. of 
pseudoscalar field leading to the relation 

(2A5A2 + A 6 (A 4 - A 3 ))(j 2 + (8A2A4 - \i)oi(J2 = -A 6 A, (62) 

whose solution is 

o 2 = Ao x + — > 0- 4 _ 2A 5 A2 + A 6 (A 4 -A 3 ) g _ A 6 A 

o\ Xq — 8A2A4 Xq — 8A2A4 

Thus in the parity breaking phase the relation between the two scalar v.e.v's is completely determined and, 
in particular, does not depend neither on p nor on p,. 

The first mass gap equation (|4"4"]) can be brought to the form 

A = 2X\of + A 5( 7ia 2 + (A 3 - Xa){oI + p 2 ) + HA{o ll p, /?), (64) 
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if one employs eq. (I5H) . Together with eq. it allows to find all v.e.v.'s of the scalar fields <Jj,p as 

functions of temperature and chemical potential. 

Let us now find the critical value of the chemical potential, namely the value where p(p c ) = 0, but 
equations (j60| . (|6T|) . (|63|) hold. Combining the two equations (|60| . ([6"Tj) 



A 6 r 2 + 4A 4 r + A 5 = 0; r = — . (65) 

CTl 

In order for a SPB phase to exist this equation has to possess real solutions. If Xq = there is only one 
solution corresponding to a second order transition, but there may exist other solutions that fall beyond the 
accuracy of our low energy model (which becomes inappropriate for small values of <7i). 

We stress that equations (|63l) and (|65l) contain only the structural constants of the potential and do not 
depend on temperature or chemical potential manifestly. Thus with the help of 

^6 

one can immediately calculate 



ot(A,\j) = \j^TA ] ^(A, A,) = r ± o£, (67) 

where of are the corresponding critical values. 

After substituting these values into equation (|64l) for each critical set of o\. 2 one derives the boundary of 
the parity breaking phase 



MA{af , Hera, Peru) = A - 2X 1 (af) 2 - X 5 afaf - (A3 - X 4 )(af) 2 . 

It must be positive at critical values of dj . The relation (|68[) defines a strip in the T — p plane where parity 
is spontaneously broken. From (|45[) one can obtain that A > and A — > 00 when T, p — > 00. It means that 
for any nontrivial solution r±, cr^, the parity breaking phase boundary exists. 

Thus we have proved that if the phenomenon of parity breaking is realized for zero temperature it will 
take place in a strip including lower chemical potentials but higher temperatures. 

B. Mass-gap equations in SPB beyond the chiral limit 

Let us now examine again the possible existence of a critical point where the strict inequality (|20p does 
not hold and for p > p cr it 



(A 3 - A 4 )cti + A 6 cncr 2 + 2A 2 [al + p 2 ^j - A 



■2 



m q d 2 



(dicri + d 2 a 2 ) F 



m q d\ 



^d 2 2 p 2 + {d 1 a 1 +d 2 a 2 ) 2 ' 
where the following consequence of equation p7| has been used: 

(7T°) diai+d 2 a 2 



F o y/d 2 p 2 + (dio-i +d 2 a 2 ) 2 
When combining equation with (TT5l) . ([T5| one finds that 



(69) 



(70) 



2 1 2 

T 2+ P 



di (X5O1 + 4A4O-10-2 + A 6 [g\ 
= 2d 2 I -A + 2X x g 2 + A 5 <7ia 2 + (A3 - A 4 ) (a 2 + p 2 ) + MA{a x ,p, /3) } , (71) 



2 1 2 

T 2+ P 



d 2 {x§(j\ + AX i cri<7 2 + A 6 [a\ 
= 2di ( - A + (A 3 - A 4 )ct? + A 6 cria 2 + 2A 2 (of + p 2 ) ) , (72) 
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where we have taken into account that o\ ^ 0. These two relations determine the v.e.v.'s of the scalar fields 
(Ti,2- If A2 = A6 = and/or A3 = A4, Xq = equations (jTTj) and (|72|) firmly fix the relation between <j\ and 
(T 2 - Otherwise an appropriate combination of these two equations still allows us to get rid of the v.e.v. of 
the pseudoscalar field [441] . Thus in the parity breaking phase the relation between the two scalar v.e.v's is 
completely determined and in particular does not depend neither on p nor on p. Using equations (jBT?)) , (ITlT) 
and (|72[) one can easily eliminate the variables p and 02, obtaining an equation for the variable <J 2 /p 2 ■ The 
latter completes the determination of the v.e.v.s. 

We notice that in the chiral limit m q — > the constants d\ , d2 become arbitrary and therefore (fTTj) , (f?2"j) 
entail three independent relations coinciding with (f6"U|), (l6"T1) , ([64"]) . 



VII. APPROACHING THE SPB PHASE TRANSITION 



Let us find the character of the phase transition to the SPB phase. For small values of p — o\ > 0, we 
know that the value of the odd parity condensate p is zero. Setting p = in equations (TT51) . (ITBI . (|4"4l and 
using (|64p and differentiating w.r.t. p we get 



k=l,2 



(73) 



or, after inversion of the matrix of second variations, 



<9 M cri = -UfcTiJp 2 -of 



V-, 



{2)a 



22 



V>ctV^ a 



< 0, 8^02 = Ufai J p 2 - erf 



(2) CT 



12 



Detl/( 2 ) CT 



(74) 



The possibility of SPB is controlled by the inequality (|20[) ; in order to approach a SPB phase transition 
when the chemical potential is increasing we have to diminish the l.h.s. of inequality (|20p and therefore we 
need to have 



(A3 — X/CjIT 2 + A6(7l(72 + 2A2CT2 



i q u, 2 



This is equivalent (using (|74)) ) to 



^A 6 cti + 4A 2 ct2 



m q d\ 



< 0. 



; )^ <(2( ^_ A4K + Am + _^_ ) 



((ilCTl + d20-2) 



(2) 



(75) 



(76) 



This last inequality is a necessary condition that has to be satisfied by the model at zero chemical potential 
for it to be potentially capable of yielding SPB. Evidently, this inequality must hold across the critical point 
in order that d^p 2 > O^d^n ) 2 > 0. 



A. Second variations in the SPB phase and character of the phase transition in the chiral limit 

Once a condensate for 7r 2 appears spontaneously the vector SU(2) symmetry is broken to U(l) and 
two charged n' mesons are expected to possess zero masses as dictated by the Goldstone theorem. For 
simplicity let us consider zero temperature. Then in the chiral limit the matrix of second variations = 



16 



(V^) ; o,6=l,2,0; reads 



l T /(2)o 



-A + 6A x o- 2 + 3A 5 o-i(7 2 + (A 3 + A 4 )o| + (A 3 - A 4 )/5 2 



/ 2 2 o 2, M + VT^l 

f-i\/ t-i — erf — 3(T 1 In - 



0"! 



1, 



(2)a 



= 3A 5 o- 2 + 4(A 3 + A 4 )cricr 2 + 3A 6 cr| + A 6 p 2 = Via, 



(2)a 
22 



1, 



-A + (A 3 + XaM + 3A 6 ctio- 2 + 6A 2 o-2 + 2A 2 p J = -V 22 , 



(2)c 



Vr 



10 

(2)<T7T 



2(1 



^4(A 3 - A 4 )(Ti + 2A 6 o- 2 )p = VioP, 
2A 6 (Ji + 8A 2 a 2 )p = V 20/ o, 



- 2 V Q0 



(2)tt 



4A 2/ o 2 = -V 00 p 2 



(2)tT 



0. 



(77) 



Now we are able to check the character of phase transition. Using consistently equations (1161) , (|44|) and the 
condition (|69p in the SPB phase one obtains the differential equations on functions <Jj(/i), p((J>), following 
the same strategy as for ([71]). 



A 2 




A 2 





V 22 Vqo — V 20 



<o, 



DetV 
VioV 2 o — Vi 2 Voo 

DetV 



DetV 



(78) 



The last derivative must be positive in order to generate parity breaking and this is guaranteed by the 
inequality (|76|) . 
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FIG. 2: The SPB phase transition of second order: the dashed line depicts the SPB breaking phase and the solid line 
stands for the v.e.v. of "dynamical" mass. The plot is only qualitative 



Let us compare the derivatives of the dynamic mass o\ across the phase transition point. For p — > ti crit — iQ 
its derivative is given by ((74)) and for ri —> p, crit + iO it is given by ([75)1 . 
Their difference reads 



duO-i 



Merit +i0 



V92V0O — Von Vn; ^ 



-»o v [ DctV DetW 2 ) CT 

(2)<r v T/ (2)<r\2 



[(Vio^'-VaoV^ 



= -4AAW M * - 22 „ ' < 0, (79) 

V 1 DetVDetW 2 )- V ; 

provided that both determinants are positive (they determine the spectrum of meson masses squared ) and 
inequality (|76[) holds across. Thus the dynamic mass derivative is discontinuous and the phase transition is 
of second order. 



B. Inclusion of current quark masses 



For non-vanishing current quark masses the deviation linear in m q in the parity breaking phase affects 
also the pseudoscalar parameters p and (n°). After usage of equation ([70)) one finds 



Acri(m g ) 

A= [ A<j 2 (m q ) ) ~2m, [ A 2 | ; (80) 
Ap(m q 

A = I Aol = ( V 



Ar\ 




A 2 




Ao 1 




P / 






/ di 


d = 






\ 4 




\ d 1 ai+d2(T2 



d\G\ + d 2 U 2 



y/d 2 2 p 2 + (d 1 a 1 +d 2 cr 2 y 



where in order to keep the leading order all parameters must be taken in the chiral limit. As to the v.e.v. 
of neutral pion field it does not need any mass corrections to the leading order and must be taken from the 
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7T°) = -arctan ( - ^ ' ) . ( si ) 



mass independent eq. p7|) 

/V°\ = — arrfian I 

v dlCTi + d 2 (J 2 

Accordingly, the mass corrections to the matrix of second variation AV^ 2 \ equation (|77|) . takes the form 

v^{m q ) = vf CT (o) + ^ a a (^ 2)ff )A a , @o = (d^d^d,)-, (82 

a=l,2,0 

C S K) = E ^ ( « 2)CT7r ) A - 

a=l,2,0 

( 2 ) 2m q di . (7T°) ( 2 ) 2m,d 2 . (tt°) 

V'o = —ft- sin TO ' = J sin ■ 



F F ' F F 

(2)tt° _ / dlOl + tfe^ (7T°) d 2 p . (it 



i\ = m,, ^ cos 7T sin ■ 



, 7 (2)7r° 2m q d 2 
V{ 2 — — ^ — cos 



Fq Fo 

1 (2W° 9, x . m„dn (7T°) 
-4 2)7r = 4A 2P 2 (0) + 16m 9 A 2 A + ^— cos ^ 

2 aiai + a 2 a 2 r 

V*r 2 ^ (7r o } 2 sin Fo ^ (7r o )2 cos (^ 



* 2 (-7T ) 

(2) ^2P . (7T°) rfl^i + d 2 Cr 2 Sm — 



^ ( 2) = 2m q dl cqs (^) 

"■2 "a 6?!^! + d 2 C2 F) 



2m 9 d 2 . (it ) 



F 
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where the r.h.s. are evaluated with the help of eqs. fp~5]) - (p~T |) . (|69| - (jTTj) and the v.e.v's for <jj,p arc taken 
in the chiral limit. We notice that convexity around this minimum implies that all diagonal elements are 
non-negative. This gives positive masses for two scalar and four pseudoscalar mesons, whereas the doublet 
of charged of tt mesons remains massless. The latter can be easily checked from the vanishing determinant 
of the last matrix V + - m ec l s - (f88]) - (|90[) . Of course, quantitatively the mass spectrum can be obtained 

only after kinetic terms are properly normalized. 

If the soft breaking of chiral symmetry occurs only in the if i channel, d 2 = then it follows from eq. ([T7| 
that light pions don't condense (ir°) and don't mix with other states as the off-diagonal matrix elements 
(|84|) . (|86|) and (l89|) vanish. The second pair of charged pseudoscalars 7r2~ becomes massless manifestly. 

VIII. KINETIC TERMS 

In this section we examine the fluctuations around the constant solutions of the mass-gap equations 
P^|) . ([TB|) . (|TT 1) . (IT%1) . and introduce appropriate notations for the fluctuations £j,n around v.e.v.'s aj, p so 
that <Zj = cfj + Sj, 7T2 = T3/0 + II. These v.e.v's <Jj,p must be used in all previous relations for the second 
variation of the potential. In calculations of the kinetic term we retain only the terms in the chiral limit 
keeping our interest to the masses of scalar and pseudoscalar mesons at the leading order of the expansion 
in current quark masses, m q . Thus in the kinetic terms we take (tti) ~ according to eqs. (IT7|) . (fT8|) . 



A. General form from chiral symmetry 

Once we have fixed the interaction to quark matter we are not free in the choice of the kinetic term for 
scalar fields. Namely one cannot rotate two fields and rescale the field Hi without changes in the chemical 
potential driver (I55[) . However the rescaling of the field H 2 is possible at the expense of an appropriate 
redefinitions of other coupling constants and this freedom can be used to fix one of the constants which 
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appear in the kinetic term. Thus we take the general kinetic term symmetric under SU(2)l x SU(2)r global 
rotations to be 



, 2 



(91) 



j,k=i 



With the chiral parameterization (|9]) one can separate the bare Goldstone boson action, 



2 2 

thin = \Yj Ajkd^d^ak + A JkV]Vktr {d^d^U} (92) 

j,k=l j,k=l 

, 2 



3 = 1 



+ ^A 22 tr {d^ 2 d^7T 2 - 2d^7T 2 ^d^rr 2 

+?d ft tfd^)(n 2 ) 2 + K+,^][7r a> ^7r 2 ]} . 

After selecting out the v.e.v. (Hj) = (oj) = <7j let us explore the kinetic part quadratic in fields. We expand 

U = 1 + i-ki/Fo H , £ = 1 + i7Ti/2F + • • • and use the notations defined at the beginning of this Section. 

Then the quadratic part looks as follows 



,(2) 



1 2 

h — An 



(93) 



-A 
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which shows manifestly the mixture between light and heavy pseudoscalar states and, in the SPB phase, 
also between scalar and pseudoscalar states. 
Let us define 



2 x 2 



(94) 



i=i 



B. Parity-symmetric phase 

In the symmetric phase p = 0, n 2 = II one diagonalizes by shifting the pion field 

< =n a - C4, (95) 
41. = ~0 /4 *«W + \{A 22 e)d^n a 2 , A 22 C 2 = > 0. (96) 

Taking into account the modification of the matrix of second variations ([85|) - ((90|) after shifting (|95"T) one finds 
the masses of light and heavy pseudoscalars to the leading order in current quark mass 



(yW*)ab = ((^) Qb -2C(^ 2) ") Qb + C 2 (^ 2) ") ab ) 

= 6 ab 2( - A + (A 3 - \,)al + A 6 a 1( x 2 + 2A 2 a 2 2 - Cm,^ 2 - + C 2 m g ^±^ ) 

= 5 ab {{A 22 -C 2 )W n + 0{ml). (97) 
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at the leading order in m q when mn 3> m w far below the P-breaking transition point in chemical potential. 
We notice that in this region the off-diagonal element does not make any influence. 

At the point of the SPB phase transition one has to impose the condition which leads to 

(V^r = S«» . U 2 - C d _l^l±Ml)\ (98) 

ai<7i + a 2 a 2 \ r Q / 

Evidently the determinant of (V^ 2 ^) vanishes and one reveals zero modes for all three pion states, one 
neutral and two charged. They represent the true Goldstone modes (in the limit of exact isospin symmetry 
m u = rrid). At the P-brcaking transition point p = 0, when taking into account the normalization of kinetic 
terms (19311 with the definitions (l94l) one finds the values of three massive modes 



2 _ A ll d 2 -2A l2 d l d 2 + A 22 d\ (dA-id) 
m u-^ detMd^+d^) - 2mq (da) ■ ( ' 

Thus in the chiral limit, at the phase transition point one reveals six zero modes and beyond the chiral limit 
only three ones. 



C. Masses of light states in SPB phase 

In the SPB phase the situation is more involved: pseudoscalar states mix with scalar ones. In particular, 
the diagonalization of kinetic terms is different for neutral and charged pions because the vector isospin 
symmetry is broken: SU(2)y — > £7(1). Namely 

In this way SPB induces mixing of both massless and heavy neutral pions with scalars. The (partially) 
diagonalized kinetic term has the following form 

42, = e^Pi* + i(i + ^-)d,i°d»i° + (A r2 - <»)a„n±3»nT 

+ i £ a^ + s^s* _ (101) 

^ . , . fn + A 22 p t + A 22 p z — ' 

j,k=l j = l 

We see that even in the massless pion sector the isospin breaking SU(2)y — > U(l) occurs: neutral pions 
become less stable with a larger decay constant. Another observation is that in the charged meson sector 
the relationship between massless it and tt' remain the same as in the symmetric phase. 

Beyond the chiral limit one can derive the masses of the lightest pseudo-goldstone states. When p m q 
then in the mass matrix ([82])- ((90)) the heavy mass parts ([82]), ([83]), jSTJ and the light mass ones ([15]), ([88]) 
- (|90[) combine into an approximately block diagonal form with additional off-diagonal elements ([84]) and 
([55]) . proportional to m q . The latter leads to factorization of the light pseudoscalar meson sector from the 
heavy meson one to the order of mi. Thus neglecting the mixture of heavy and light states one deals with 
the light sector built of ([85]) . (|88]) - ([90]) which after diagonalizing the kinetic term by (|100[) projected on the 
light state sector gives the light pseudoscalar masses 



„ , A 2 iP 2 \- x fd 1 a 1 +d 2 a 2 (tt ) d 2 p . 
TO#o =2m^l+-^J ( - 2 cos— --jjg-ffln- 



F 2 F F 2 F 



mi ± = 0, (102) 
2 _ cosi fr {, ,d 1 a 1 + d 2 a 2l _ (tt ) 



TOtt± = 2m a —. 7777-— : d 2 — C tan ^ • (103) 

n 9 (A 22 -C 2 )(rfiCTi+rf 2 a 2 ) V (t°> j 



Thus in the SPB one finds two massless charged pseudoscalars and three light pseudoscalars with masses 
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FIG. 3: Masses of light states in SPB phase: the vacuum masses are taken from pcj /. The light and heavy charged 
pseudoscalars are depicted with solid lines, the dotted line corresponds to the neutral light pseudoscalar and the 
dashed line stands for the neutral heavy one. The plot is only qualitative. 

linear in the current quark mass. These equations represent the generalization of the Gell-Mann-Oakes- 
Rcnner relation in the phase with broken parity. 

We notice that the masses of neutral and charged pseudoscalars do not coincide in the well developed SPB 
phase, just realizing the spontaneous breaking of isospin symmetry. One can also guess that the manifest 
breaking of SU (2) symmetry due to different masses of u and d quarks will supply the Goldstone bosons tt^ 
with tiny masses proportional to the difference m u — m^. 

IX. DESCRIPTION OF NUCLEAR MATTER 

The normal density of infinite nuclear matter [26| is go ~ 0.154-0.16 fm~ 3 that corresponds to the average 
distance 1.8 -j- 1.9 fm between nucleons in nuclear matter. We assume that the quark matter is equivalent 
to nuclear matter when their average baryon densities coincide, at least in what respects meson properties. 
This normalization we will apply at the normal baryon density while noting that there also exist different 
proposals [4l| to confront quark and nuclear matter for higher densities. 

Thermodynamical characteristics are the pressure, P, and the energy density, e. The pressure is determined 
in the presence of chemical potential by (l49t , defined for aj satisfying the mass gap equation. Evidently the 
pressure at zero nuclear density vanishes. In this case the energy and baryon densities are related to the 
pressure as follows 

e = -P + N c hq b ; d^P = N cQB - d e e = N ct i. (104) 
The direct connection between energy density and pressure reads 

P = Q%d e (j^j ■ (105) 

Evidently the energy per baryon has an extremum when the pressure vanishes. Since the pressure is an 
increasing function of the density as we have seen, obviously vanishing at zero density, and infinite nuclear 
matter is stable (thus implying zero pressure) the phase diagram in the P,Qb plane is necessarily discon- 
tinuous with values of density in the interval (0, go) not corresponding to equilibrium states (go is nuclear 
matter density). We will see below how this is realized in our model. 

A. On the way to stable nuclear matter 

Our model consisting of two scalar isomultiplets is still somewhat too simple in one aspect. The stabi- 
lization of nuclear matter requires not only attractive scalar forces (scalars) but also repulsive ones (vector- 
mediated). Conventionally (26|, the latter ones are associated to the interactions mediated by the iso-singlet 
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vector uj meson. Let us supplement our action with the free uj meson lagrangian and its coupling to quarks 



AC 



1 



1 



(106) 



with a coupling constant g^qq ~ 0(l/y/N c ). After bosonization of QCD, on symmetry grounds, any vector 
field interacts with scalars in the form of commutator and therefore lo^ does not show up in the effective 
potential Hj fields to the lowest order. However in the quark lagrangian the time component ujq interplays 
with the chemical potential and it is of importance to describe the dense nuclear matter properties. Let 
us assign a constant v.e.v. for this component ^^(wo) = Co, Then one needs to compute the modification 
of the effective potential due to the replacement p — > ji + Co = p. The variable Co, and accordingly p, is 
dynamical and in addition appears quadratically in the mass term in (|106[) which reads 



1 



1 (A - p) 2 

2 Gui 



Q 2 - 1 

iluqq ^ (n(__ 

ml ~ { N r _ 



(107) 



The term (I107[) supplements the effective potential ([56)) : V e g (p) = V e g (p) + AK;(/i,/i). Correspondingly 
the extremum condition for the variation of the variable p involves both the scalar part of the effective 
potential (|5"6"|) and the vector one (|107|) and due to ([5"2"j) takes the following form 



N c Qb{p) 



(a) 



(108) 



from this one finds p(p) after solving the mass-gap equations (J3U), (jTHJ) and (ITBl . 
Finally the extended effective potential at a minimum reads 



V edri = -2 A ( <J i(/ i )+ (7 2(A) 



V(A) 



/ \3/2 8_/y2 , \ 3 i / \ 

(A 2 -^i(A)J + G u —[p 2 -(7i (A) J J0(a-<ti(A)J 



(109) 



Let us define the v.e.v. 's of scalar field o\ in vacuum at the two minima as er*(0) < crj(0) . Let us select 
out the parameter subspace such that the minimum corresponding to a\ (0),cr|(0) is lower than the another 
minimum at cr* (Q) , a\ (0) . Then for parity-even matter p = 0, one seeks for the nuclear matter stability 
at a value of chemical potential p s with cr^(0) < p s < crj(0) . The corresponding baryon matter stability 
condition AP = P(a{(0)) - P(af(p s )) = 0, eq. (j49]), can be formulated as 



g^nIqI^s 



(110) 



taking into account (| 109[) and (|107[) . Herein /} s is related to the physical value of p s by (llOSp and it is 
assumed that parity is not violated p = 0. This relation represents the condition for the formation of stable 
symmetric nuclear matter in result of first-order phase transition [26| . It can be fulfilled by an appropriate 
choice of the vector coupling constant G^ as typically the first term in the r.h.s. of (| 1 10[) is smaller than the 
one on the l.h.s. 

At finite temperatures one has to modify the thermodynamic relations. The modification of the effective 
potential due to to mesons is given by (|107[) . Thus V e g (p,(3) = V e Q (p,(3) + AV L0 (p,p) which should 
henceforth be used in all the previous thermodynamical formulae. The replacement p p makes all 
expectation values depend rather on p which is determined via the variation of V e Q 



Gu 



= -N cQB (p,p,cj 1 ) =dnAV eS (fi,p) 



(111) 



The saturation point at p = p s where nuclear matter forms is characterized by the energy crossing 
condition for P, T f 0, 



+\t(s[P, - e(fi, p s ,4)) + GcjN 2 (g 2 B {(3, p s ,a* x ) - g%(p, p s , tr})) , 



(112) 
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where fi s is related to the physical value of fj, s by equation and = cr(/i s ,/3); a* = cr*(fl s ,f3). This 

relation represents the condition for the existence of symmetric nuclear matter. It can be always fulfilled by 
an appropriate choice of G u . It can be shown that normal nuclear matter is formed at the chemical potential 
/i s ~ 303MeV stabilized by u> meson condensate with G^ ~ (10 -j- 15)GeV~ 2 in a satisfactory agreement to 
what is known from other model estimations [421. 



B. First-order transition without chiral collapse 



A viable model of dense baryon (quark/nucleon) matter must reveal the phase transition to a stable 
bound state at the normal nuclear density qb = go for infinite homogeneous symmetric nuclear matter at 
the so called "saturation point" . This phase transition is believed to be of first order similar to the vapor 
condensation into liquid: from droplets of heavy nuclei to a homogeneous nuclear liquid. However in simple 
quark models of the NJL type |38[ this phase transition (for vanishing current quark masses) goes to the 
chirally symmetric phase with zero dynamical mass (zero v.e.v. of scalar fields), so called "chiral collapse". 
When it happens the typical baryon density is substantially larger that the normal one qb, c — 2.8go- For this 
reason these simple models cannot be a reasonable guide to phase transitions in very dense nuclear matter. 

Let us examine under which conditions the saturation point in our model happens to be at normal nuclear 
density and is not accompanied by the chiral collapse keeping the dynamical mass different from zero. To 
analyze this problem we have to examine the pressure in cold (T = 0), dense baryon matter. 

We remind that in so far as our system undergoes spontaneous CSB the effective potential (|3"2")l does 
not reveal any minimum at the origin in variables Oj (for \i = 0) and may have either a saddle point, 



det 



y( 2 ) 



((Tj = 0) < or a maximum, det 



(0) > 0, tr 



< 0. One has a positive def- 



inite matrix of second variations (|2~Tj) . (|22|) of the effective potential in the vicinity of a CSB solution, 



det 



t>( 2 ) 



(oj) > 0; tr |l>( 2 )| (crj) > 0. It means that V eS (aj) < V eS {0). 



These properties allow us to guess that at some value of chemical potential fi s < a\ and smaller values 
of v.e.v. 's for scalar fields a\ < fi s the deficit in scalar background energies on the left-hand side of f|l 10[) 
may be exactly compensated by contributions from the nuclear density and omega-meson repulsion on the 
right-hand side so that -P(c*, (i s ) = and the system undergoes a first-order phase transition to the stable 
quark (nuclear) matter. 

Let us now prove that for a large variety of coupling constants admitting CSB (and SPB, see below) one 
of the v.e.v. a* ^ in the chiral limit, and the chiral collapse is impossible. Indeed, suppose that a* = at 

/i* < o\ where the pressure vanishes then 



3 



v eS (4/i = o) 



In this case the second variation matrix for the effective potential (f21jl . (|58|) reads 



-Vtf )a = -A 
2 11 



-A/V) 5 



Vg>° = 0, \v™> = -A. 



(113) 



(114) 



In order to induce SPB one takes A > (see below). Then from (|114p one finds that for any value of /i* 
the second variation matrix is never positive definite and one reveals either a saddle point or a maximum at 
a presumed saturation point whereas a maximum remains for vacuum values with fj,* < a\ . As we have to 
guarantee the existence of stable nuclear matter with normal baryon density we consider further on A > 0. 

It is instructive to reduce the two-multiplet sigma model to a one-multiplet lagrangian associated to a NJL 
quark model. In relations (|114[) it simply corresponds to taking only one matrix element V^ a to describe 
the behavior around extremum at the origin. Evidently, there is always a value of \i* for which it becomes 
positive and chiral collapse is inevitable. But when comparing with our model one concludes that the reason 
for appearance of chiral collapse is not the absence of confinement [38| but the inefficiency of a one-channel 
linear sigma model in representing the complicated chiral dynamics in hadronic physics. 



C. (In)compressibilities: matching quark and nuclear matters 



The (in)compressibility in the quark matter must be defined as K(p) = d gB P, where the derivative is 
made with the help of the function Qb{^) given in ([55]) . For a zero pressure state (such as stable nuclear 
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matter) this equals 



K = S 2 B dl(j-)\ p=Q , fUo, 



which must be positive since it corresponds to a minimum of the energy per baryon. In our model this is 
indeed the case 

K(u) = N cQB d eB fx = r + 9G u g B (n), d m = -AMaJp - a\ -^ij- < 0. (116) 



where N c — 3 is assumed. 

In order to adjust the properties of stable baryon matter in a quark description we parameterize the 
(in) compressibility in the following way: K — aN c g B d g fi . Let us show that the matching of quark 

and nuclear matter at the saturation point is provided by the normalization factor a = 9 in meson-nucleon 
models and a = 1 for quark-meson models. The incompressibility must be positive (giving a minimum of the 
energy per baryon) when a stable nuclear matter is formed at zero pressure. The derivative of the dynamical 
mass is given in (|116j) wherefrom it follows that the derivative of the pressure is always positive. Thus if 
there is a solution with a\ = ci(/i*) < \i* < a\ providing P = 0, the phase transition emerges to the stable 
nuclear matter state. 

In the terms of nuclear d.o.f. one defines 

£=-P + ^NQB- (117) 

At the saturation point P = and 

^ijv = — = energy per baryon = itln — Eb oun d — (939 — 16)MeV — 923MeV. (H8) 

QB 

The quark matter chemical potential is defined as d e e = /j-n — N c fi. Therefore at the saturation point 
fx s = 30&MeV. Let us use this point for the quark-hadron matching 

Qb(vn) = Qb(^s), Pf.n = (^n) 2 - {m N ) 2 ~ p 2 Fq = (us) 2 - (cr*) 2 \ ( 119 ) 

if we neglect the vector meson shift fi ~ /i. However the matching of densities does not provide the equivalence 
of derivatives w.r.t. chemical potentials. Namely around the saturation point 

|^ ~ 2fi N = 2N ct i s 2 Ms , (120) 

where (subdominant) derivatives of dynamical masses are neglected for a moment. Now let us try to extend 
the matching to the (in)compressibilities, providing the correct N c factors. For hadron matter 

\OUnJ UN NrUs 



whereas for quark matter, 



they match each other for N c = 3. One could do things even better. If the coefficient for hadron matter we 
redefine 9 — > 3N C and we introduce the coefficient 3/7V c for quark matter then both definitions match for 
any iV c to the leading order. At least one then could succeed in their matching around normal density. 
Going back to quark matter description it would mean that 

K((m s ) = + 3N c G u g B (» s ), (123) 

has a finite limit at large A^ c coinciding with what we get for hadron matter if one remembers that G u ~ 1 /N c 
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D. Saturation point meets spontaneous parity breaking 



Let us search for the domain of parameters in the model providing the realization of both stable nuclear 
matter and the regime of SPB. The former is associated with a first-order phase transition and implies the 
existence of two minima at zero chemical potential which are possibly moving when the chemical potential 
increases. The highest, metastable minimum must start moving at chemical potentials fi smaller than the 
value of the dynamical mass of the lowest minimum, a\ and larger than the v.e.v. er* at the highest minimum, 
a \ < A* < a \ ■ This metastable minimum may reach the lowest one if the density and omega meson effects 
are taken into account. Then a first-order phase transition to normal nuclear matter occurs when pressures 
become equal, eq. (|110|) . 

In order to simplify our search we make a particular choice of A5 = (not a reduction by (|29|) 0. In this 
case one of the solutions is = and 2Ai(<7 1 °' ) ) 2 = A and it is a minimum as it follows from ([2~Tj) . (l22l 

(in the chiral limit) provided that Ai > 0, (A3 + A4) > 2A X . When — a higher symmetry Z 2 X Z 2 
arises for the effective potential in the vicinity of such a minimum as the contribution of the vertex with Ag 
into the second variation vanishes with ■ For 02 7^ one can obtain eq. (|35|) for the ratio t = u-il^x- As 
it is analyzed in subsection 4.2, it has, in general, one or three real roots. For our purposes eq. (|35|) must 
have three real solutions: one corresponding to a minimum > and two corresponding to saddle points 
< 0,t^ > 0. The inequality controlling the existence of three real solutions is derived in subsection 4.2 
from the analysis of the minimum of the polynomial (1351) . Finally for a given solution t' 3 ' one finds a unique 
solution for <j\ > from JT5) (or (JUJ) and {H 



Let us assume the minimum with a 2 — to be the higher one at zero chemical potentials = a\ . For 

this choice to be realized it is sufficient to fulfill the inequality <j\ = u\ < a\ = a\ . It turns out that in 
order to provide it one has to satisfy the inequality 

2A 2 t 2 + ^A 6 t+(A 3 + A 4 -2A 1 ) < 0, (124) 



which implies 



9A 2 ,>32A 2 (A 3 + A 4 -2Ai); < i (3) < -— -. (125) 

4A2 



When at the critical value of fi = /i s < ctJ, the solution with a\ = describes a saturation point then the 
further evolution of the meson background for higher chemical potentials is characterized by the following 
equation a*(/j,*) 

2Ai(cti) 2 (/i) = A-J\fA(al,fjL). (126) 

As the last term is monotonously increasing with chemical potential the v.e.v of scalar field is decreasing. 
Now we approach to the P-breaking regime and employ eq. (|S(J]) at the expected phase transition point. Its 
solution is 

A3 — A4 l\\ 



Thus the feasibility of spontaneous P-breaking depends on the realization of the inequality o\ iC < a\ < a\ 



# 




FIG. 4: Saturation point meets spontaneous parity breaking: at li — fi s the pressures for the two solutions a,a* 
become equal and the solutions interchange realizing the 1st order phase transition. At a larger chemical potential 
fx c the 2nd order SPB phase transition occurs. 



Let us collect the inequalities providing the required convexity of the two minima and the very existence 
of both the stable nuclear matter and a parity breaking phase for higher densities (see [TH, [l9j]) 



3Ag 
~4A 2 



> > max 



3Ae. 
8A 2 ' 



4A 4 



^1,2,3,4 > 0, A 3 > A 4 , -A 6 > 8A 2 A 4 > A 



A > 0, (A 3 ± A 4 ) > 2Ai, (A 3 + A 4 ) > 2A 2 



(128) 



in addition to those ones derived above. For a more definite numerical estimation of these six constants there 
is not at present enough experimental or phenomenological information, although it can be shown that the 
tentative values assumed in [18| for Ai ~ 0.15, A3 ~ 4, A ~ 0.03Gel^~ 2 may lead to the occurrence of SPB 
at about three times normal nuclear densities. 



X. CONCLUSIONS 



Let us summarize here our main findings. Parity breaking seems to be quite a realistic possibility in 
nuclear matter at moderate densities. We have arrived at this conclusion by using an effective lagrangian 
for low-energy QCD that retains the two lowest lying states in the scalar and pseudoscalar sectors. We 
include a chemical potential for the quarks that corresponds to a finite density of baryons and investigate 
the pattern of symmetry breaking in its presence. We have found the necessary and sufficient conditions for 
a phase where parity is spontaneously broken to exist. In general this phase is bound and it extends across 
a range of chemical potentials that correspond to nuclear densities where more exotic phenomena such as 
color-flavor locking or color superconductivity may occur. It also extends to finite temperature. 

Salient characteristics of the phase with parity breaking would be the spontaneous breaking of the vector 
isospin symmetry SU(2)y down to U(l) and the generation two additional massless charged pseudoscalar 
mesons. We also find a strong mixing between scalar and pseudoscalar states that translate spontaneous 
parity breaking into meson decays. The mass eigenstates will decay both in odd and even number of pions 
simultaneously. Isospin breaking can also be visible in decay constants. 

Let us mention here several possible experimental or phenomenological signatures of parity breaking, 
discussed previously in [l8|, [l9[ in the chiral limit. 
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a) Decays of higher-mass meson resonances (radial excitations) into pions. Resonances do not have a 
definite parity and therefore the same resonance can decay both in two and three pions (in general into even 
and odd number of pions). 

b) At the very point of the phase transition leading to parity breaking one has six massless pion-like states 
in the chiral limit and the two massless charged states when the quark masses are taken into account. After 
crossing the phase transition, in the parity broken phase, the massless charged pseudoscalar states remain 
as Goldstone bosons enhancing charged pion production, whereas the additional neutral pseudoscalar state 
becomes massive. 

c) Reinforcement of long-range correlations in the pseudoscalar channel. 

d) Additional isospin breaking effects in the pion decay constant and substantial modification of F v ' for 
massless charged pions, giving an enhancement of electroweak decays. 

We think that our conclusions are drawn in a region of parameters where our effective lagrangian is 
applicable and, while obviously we cannot claim high accuracy in our predictions, we are confident that the 
existence of this novel phase is not an spurious consequence of our approach but a rather robust prediction. 
It would surely be interesting to investigate how this new phenomenon could possibly influence the equation 
of state of neutron stars (the density of such objects seems to be about right for it). 

One could hope that lattice methods could shed some light on this issue and confirm or falsify the existence 
of this interesting phase in dense nuclear matter. The problem, as discussed in the introduction, is far from 
trivial. One could probably use the expansions for small values of (i at finite t emp eratures to check some of 
our expressions. For this matching the natural approach is a hot hadron gas [43| . Conversely, it would be 
possible to extend our techniques to the case of isospin chemical potential, that is physically uninteresting 
but amenable to numerical simulations as the fermion determinant is positive in this case (of course there is 
no true parity violating phase in this case[l2j]). 

We have seen that SPB appears to be a generic and rather robust possibility offered by QCD at finite 
density and temperature. This conclusion has been reached by making sure that the low-energy modelling 
of QCD is as accurate as it is reasonable to expect. We have seen that as a bonus we get a rather good 
description of several aspects of nuclear physics; in particular a good description of the physics associated to 
the condensation transition where nuclear matter becomes the preferred solution. The model is rich enough 
to provide the relevant characteristics while avoiding some undesirable properties of simpler models, such as 
the chiral collapse. Other nuclear properties such as (in) compressibilities are well described too. 

We have presented our results trying to avoid as much as possible specific numerical values for the different 
quantities and parameters. Not only is this procedure more general but also the logical connections are better 
outlined. The known restrictions and meson properties of QCD are usually included as inequalities. The 
main conclusion of our studies is that spontaneous parity breaking may be a rather generic phenomenon at 
finite density. 
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